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Abstract 

We propose to use Fermi-Dirac distributions for quark and antiquark 
partons. It allows a fair description of the x-dependence of the very recent 
NMC data on the proton and neutron structure functions Fr?(x) and F£(x) 
at Q 2 = 4 GeV 2 , as well as the CCFR antiquark distribution xq(x). We show 
that one can also use a corresponding Bose-Einstein expression to describe 
consistently the gluon distribution. The Pauli exclusion principle, which 
has been identified to explain the flavor asymmetry of the light-quark sea 
of the proton, is advocated to guide us for making a simple construction of 
the polarized parton distributions. We predict the spin dependent structure 
functions gf(x) and g™(x) in good agreement with EMC and SLAC data. 
The quark distributions involve some parameters whose values support well 
the hypothesis that the violation of the quark parton model sum rules is a 
consequence of the Pauli principle. 
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1 Introduction 



Many years ago Feynman and Field made the conjecture' 1 ' that the quark 
sea in the proton may not be flavor symmetric, more precisely d > u, as 
a consequence of Pauli principle which favors dd pairs with respect to uu 
pairs because of the presence of two valence u quarks and only one valence 
d quark in the proton. This idea was confirmed by the results of the NMC 
experiment' 2 ' on the measurement of proton and neutron unpolarized struc- 
ture functions, F 2 (x). It yields a fair evidence for a defect in the Gottfried 
sum rule' 3 ' and one finds 



J ^W(ar) - F2{x)} = 0.240 ± 0.016 (1) 



- 1 dx , 
/o x 

instead of the value 1/3 predicted with a flavor symmetric sea, since we have 
in fact 

I G = \{u + u-d-d) = - + -(u-d). (2) 

A crucial role of Pauli principle may also be advocated to explain the well 
known dominance of u over d quarks at high x,' 4 ' which explains the rapid 
decrease of the ratio F£(x) / 'F^ix) in this region. Let us denote by q*(q^), u 
or d quarks with helicity parallel (antiparallel) to the proton helicity. The 
double helicity asymmetry measured in polarized muon (electron) - polarized 
proton deep inelastic scattering allows the determination of the quantity 
Ai(x) which increases towards one for high a;,' 5 ' suggesting that in this region 
w dominates over u^, a fortiori dominates over d) and o^, and we will see now, 
how it is possible to make these considerations more quantitative. Indeed at 
Q 2 = the first moments of the valence quarks are related to the values of 
the axial couplings 



t i , t l+F- D , l-F + D . , 
u\ al = l + F, v* nl = l-F, < al = , < al = , (3) 

so by taking F — 1/2 and D = 3/4 (very near to the quoted values' 6 ' 0.477 ± 
0.011 and 0.755 ± 0.011) one has uL = 3/2 and uL = 1/2 which is at the 
center of the rather narrow range (c£, 4ai) = (3/8^/8). The abundance of 
each of these four valence quark species, denoted by p v3j \, is given by eq. (Q) 
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and we assume that the distributions at high Q 2 "keep a memory" of the 
properties of the valence quarks, which is reasonable since for x > .2 the sea 
is very small. So we may write for the parton distributions 



p{x) = F(x,p val ) (4) 

where F is an increasing function of p V ai- The fact that the dominant distri- 
bution at high x is just the one corresponding to the highest value of p va i, 
gives the correlation abundance - shape suggested by Pauli principle, so we 
expect broader shapes for more abundant partons. If F(x,p va \) is a smooth 
function of p va \, its value at the center of a narrow range is given, to a good 
approximation, by half the sum of the values at the extrema, which then 
implies^ 

u^{x) = -d(x). (5) 

This leads to 



Au(x) = uHx] 



uHx) 



u(x) — d(x) 



(6) 



which allows to relate the contribution of u quark to gi(x), the proton po- 
larized structure function, to the contributions of u and d quarks to F%(x) — 
F?(x) i.e. 



X 9l( X y 



(F$(x) - F?(x)) 



u+d 



Moreover, since Ad v 



■ F-D 



Au v 



— l/AAu v and e\ 



small negative correction to Qi(x) from d quarks and therefore 



(7) 

l/4e^, we expect a 



xgl{x) 



(F$(x) 



(8) 



at least in the x region dominated by valence quarks. This is confirmed by 
experiment ' 2 ' 5 ^, which shows good agreement with eq. for x > .2, with a 
slight overestimate of gx(x) at small x, partly due to the fact that we have 
neglected the d quark contribution Ad, which is expected to be negative ™. 
All these features are encouraging to take rather seriously the importance of 
Pauli principle and have led us to propose a new description of the quark 
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parton distributions in terms of Fermi-Dirac expressions as functions of the 
scaling variable x according to' 8 ' 

x 

Here x(p) plays the role of the "thermodynamical potential" for the fermionic 
parton p of a given flavor and helicity, x is the "temperature" and f(x) is 
a weight function whose support is the open range (0,1). As we will see, 
we will be able to reproduce the NMC data' 9 ' for F$(x) and F£(x) taken at 
Q 2 = 4 GeV 2 , and the antiquark data from neutrino deep inelastic scattering 
P-°l. This will allow the determination of the small number of parameters 
entering in eq. (|). Keeping the same value for x, we will also produce a 
gluon distribution G(x) in terms of a Bose-Einstein expression similar to 
eq. (^) and consistent with what is known from experiment' 11 ' 12 '. We will 
find for x(u^) a much larger value than for x(u^) (i.e. x(u^) ~ 4x(u^)) which 
suggests that the defect in the Gottfried sum rule is a consequence of the 
fact that the u u* bus" is almost full and therefore the sea quarks have to fill 
up all available space in the "w^ bus" and in the u d bus" . This allows us to 
propose an interesting solution of the spin crisis, where we will have a fair 
description of g\{x) obtained by EMC' 5 ' and also of the neutron polarized 
structure function g™(x) obtained very recently at SLAC' 13 '. 

In the next section we present our new approach for parton distributions 
and we will specify our parameterization for u and d quarks, for antiquarks 
and for gluons. In section 3, by using the most recent and accurate avail- 
able low Q 2 unpolarized deep inelastic data, we will see that it is possible to 
complete the determination of all these distributions in terms of eight free 
parameters and we will also present our predictions for polarized deep inelas- 
tic scattering. In section 4, we will discuss shortly the outcoming results of 
this approach for some quark parton sum rules. We will give our concluding 
remarks in section 5. 

2 New Approach to Parton Distributions 

Let us first consider u quarks and antiquarks and we assume that m t (x), 
u^(x), vJ(x) and u^(x) are expressed in terms of Fermi-Dirac distributions of 
the form given by eq. (|9|). The thermodynamical potential x(p) is a constant 
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parameter which is expected to be different for each parton p, whereas the 
temperature x can be a universal constant and we will take for fix) a simple 
formf] 

f{x) = Ax a . (10) 

Concerning the d quarks and antiquarks, following our above arguments we 
will assume that 

d(x) = ^ (5') 

which is a slight modification of eq. (H), to account for the fact that F 
is not exactly 1/2. As also indicated above, Ad(x) is not very large and 
moreover the potentials associated to d^(x) and d^(x) are expected to satisfy 
the following constraints 

< x(d^) — x{u^) ~ a;(ir) — x(dJ) < x. (11) 

Therefore this observation justifies, to a reasonable approximation, the fol- 
lowing choice 

M{x) = ~ kfl - x \ : ,-ZT tl y (12) 

an expression with no new parameter, except the normalization factor k 
which will be fixed by requiring for the first moment 

Ad = Ad val = F-D. (13) 

We therefore assume that the d sea quarks are not polarized^ and consistently 
we also take for the d antiquarks, 

d\x)=d i (x)=u i {x). (14) 



*With this choice p(x) does not go to zero when x — ► 1, but it has a fast decrease 
coming from the exponential in the denominator. 

t Clearly whereas one needs a large negative polarization for u sea quarks and anti- 
quarks, it is less crucial for the case of the d's which where assumed to be slightly polarized 
in refs.[14] and [15]. 
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Concerning the strange quarks we first take in accordance with the data' 10 ] 



/ x — / \ u(x) + d(x) 
s(x) = s{x) = . (15) 

Finally for the gluon distribution, for the sake of consistency, one should 
assume a Bose-Einstein expression given as 

r>t \ 16 f( x ) / 1fi \ 

X 

with the same temperature x, a specific potential x(G) and where the factor 
16/3 is just twice the ratio of the color degeneracies of gluon and quarks since 
G(x) is the unpolarized gluon distribution. 

This completes the parameterization of all parton distributions we will 
use and we will now proceed to the determination of the eight parameters 
we have introduced by means of the description of low Q 2 unpolarized deep 
inelastic scattering data. We will also give the predictions it leads to for the 
polarized structure functions g{(x) and gi(x). 



3 Unpolarized and Polarized Deep Inelastic 
Scattering 

To determine our parameters we have used the most recent NMC data^Q on 
F$(x) and F^x) at Q 2 = 4 GeV 2 , together with the most accurate CCFR 
data on the antiquark distribution [10]. 
We get the following set of parameters 

x = 0.132 , A = 0.579 , a = -0.845, (17) 

x(u r ) = 0.524 , x(u L ) = 0.143 , £(m t ) = -0.216 , x(u l ) = -0.141. 

For the fraction of the total momentum carried by quarks and antiquarks we 
find 



tin ref.[9] the result quoted in eq.(l) has been slightly changed to Iq = 0.258 ± 0.017. 
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/ xu{x)dx = 0.278 , / xd{x)dx = 0.075, 
Jo Jo 

f 1 x\u{x) + d(x) + s(x) + s(x)]dx = 0.084. (18) 
Jo 

We show the results of the fit for F%{x) — F£(x) and F^x)/ F% (x) in Figs. 1 
and 2 and in Fig. 3 our prediction for the distribution of the antiquarks com- 
pared to the data' 10 ]. The behavior of the antiquark distribution is correctly 
reproduced for x > 0.1 but for small x values we don't have the fast increase 
shown by the data. This is due to our simplifying assumptions that the x's 
and x are taken to be constant. Actually when x — > 0, from Pomeron uni- 
versality, one expects xu(x) = xd(x) ^ 0, contrary to the present situation 
where we took for quarks and antiquarks a universal f(x) such that xf(x) 
vanishes when x goes to zero. Clearly this is not adequate for antiquarks 
and could be improved in a more sophisticated version of this approach for 
example by allowing the x's to depend on x. In particular, we underestimate 
u and d so we don't satisfy very accurately the obvious constraints 

u — u = 2 and d — d=l. (19) 

We now turn to the gluon distribution (see eq. (Jf6|) whose only free pa- 
rameter x(G) has been fixed by the requirement that gluons carry the fraction 
of the total momentum not carried by quarks and antiquarks, that is 0.563, 
according to eq. flilf) . We find a very small value x{G) = —0.012 and we 
display in Fig. 4 our prediction which is fairly consistent with some prelim- 
inary indirect experimental determination from direct photon production' 11 ] 
and from neutrino deep inelastic scattering' 12 ' at Q 2 = 5 GeV 2 . 

Going back to quark distributions we can now test our approach by look- 
ing at the predictions we obtain for the polarized proton and neutron struc- 
ture functions. So far, all the parameters have been fixed, except the nor- 
malization factor k of Ad (see eq. (0)) which was found to be k = 0.787 
after imposing eq. (|j"3"|). The results at Q 2 = 4 GeV 2 are shown in Fig. 5 for 
g\{x) and in Fig. 6 for gi(x). For g\{x) the agreement with the data which 
corresponds to (Q 2 ) = 10 GeV 2 , is very satisfactory but no significant Q 2 
dependence has been observed by EMC [5] . For g™(x) we have the general 
trend of the data at (Q 2 ) = 2 GeV 2 , which is small, because there is a strong 
cancellation between Au(x) weighted by 1/9 and Ad(x) weighted by 4/9. 
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However we fail to reproduce accurately the structure around x = 0.1, in 
spite of the fact that the driving contribution Ad(x) has a minimum at the 
correct value because x(u^) = 0.143 (see eq. (|T2"|)). These predictions are 
very encouraging and should certainly be improved. 



4 Quark Parton Model Sum Rules 

There is an intriguing consequence of the hypothesis that the defect in the 
Gottfried sum rule follows only from Pauli principle. In this framework since 
in the "v} bus", there is less available space than for the other partons, we 
expect that the defect in the Gottfried sum rule arises from the fact that the 
sea must contribute less to the u* distribution. This implies from eq. (|14"D 

An = Au val + u-d = 2F + u-d = 0.814 ± 0.046, (20) 

which is smaller than Aii va i because d > u. Let us now consider the Bjorken 
sum rule' 16 ^ which reads 



ri 1 
/ _-J n = / \gf(x) -g?(x)]dx= -(Au + Au- Ad- Ad) 
Jo 6 

= \{F + D) = ^(A« val -Ad val ) (21) 
6 6 

where we have used eq. (|3]) for the last equality. Since in our approach 
Ad + Ad = Adyax = F — D, in order to satisfy the Bjorken sum rule we 
should have Am + Am = Aw va i = 2F and therefore Au should be positive at 
variance with Aw sea . Of course, given the fact that As = As = 0,' 17 l this 
leads for the parton Ellis- Jaffe sum rule' 18 ' 



I P = I 9i{x)dx = -(Au + Au) + — (Ad + Ad) 
Jo 9 18 

f n 

0.196 ±0.006 (22) 



2 U 

to be compared with the EMC result 0.126 ± 0.015 ± 0.009. However if one 
chooses 
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Am = 2F + u - d, Au = Aw sea = u - d, Ad = F-D, Ad = 0, (23) 



one gets 



and 



FDA, - 
I p = -- — + -(u-d) = 0.134 ± 0.015 
2 lo y 



I n = \f — ^-D + -(u - d) = -0.024 ± 0.009 (24) 
3 9 9 

in fair agreement with the proton EMC result ^ and with the neutron SLAC 
result^ which is I n = -0.022 ± 0.011.0 

This choice (eq. which is obtained by assuming the crucial role of 

Pauli principle for parton distributions implies that the theoretical framework 
for the derivation of the Bjorken sum rule should be reconsidered. 



5 Concluding remarks 

We have obtained a good description of the main features of unpolarized 
and polarized deep inelastic scattering by using Fermi-Dirac distributions for 
quark and antiquark partons. This supports the idea, previously proposed^ , 
that Pauli exclusion principle plays a crucial role in deep inelastic phenom- 
ena and leads to the intringuing conclusion that partons interact incoherently 
but are bound to obey Fermi-Dirac statistics in the x variable. Concerning 
the parameters of these distributions, we have found a universal temper- 
ature x for quarks, antiquarks and gluons, whereas the thermodynamical 
potentials x(p) are positive for quarks, negative for antiquarks and small for 
gluons. As expected, we find that x(u*) dominates over all the other ones, 
so the " bus" is rather full which is naturaly related to the defect of some 
quark parton model sum rules. In its present form, our approach fails to 
reproduce the fast increase of the antiquark distribution at very low x and 
this deserves some further comments. One obvious reason is that within 

§A much less accurate CERN experiment,! 19 ! leads to a different result, i.e. /„ = 
-0.08 ±0.04 ±0.04. 
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the parton model approximation the distributions are depending on x only, 
not on Q 2 , and x and x(p) were taken to be constants. Instead, we could 
have assumed that they depend also on the final hadronic invariant mass 
(P + q) 2 = M 2 + Q 2 ((l — x)/x) which varies rapidly at small x. We may 
observe that this fast increase of the distributions at small x was advocated 
in the framework of the multiperipheral approach of deep inelastic scatter- 
ing, where a singular behaviour x~ 3 ^ 2 was suggested' 20 ' and this trend was 
recently observed in the data' 21 ' for the first time. In this case one can speak 
pictorially of a liquid of partons for which the Fermi-Dirac gas description is 
inadequate. The parton distributions contain two phases, a gas contributing 
to the non singlet part which dominates at moderate and large x and a liquid 
contributing to the singlet part which prevails at low x. This is not a new 
situation and it corresponds to the well known two components picture, with 
ordinary Regge trajectories related to resonances and the Pomeron. Finally, 
we would like to add that we have been recently aware of a work^ 22 ] where 
deep inelastic structure functions are calculated in a statistical model of the 
nucleon considered as a gas of quarks and gluons in the framework of the MIT 
bag model. Despite the commun reference to quantum statistics, the two ap- 
proaches are rather different and in ref. [22] only unpolarized distributions 
have been considered. 
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Figure Captions 



Fig.l The difference F% (x) — F^x) at Q 2 = 4 GeV 2 versus x. Data are from 
ref.[9] and the solid line is the result of our fit. 

Fig.2 The ratio F^{x)/F^(x) at Q 2 = 4 GeV 2 versus x. Data are from ref.[9] 
and the solid line is the result of our fit. 

Fig. 3 The antiquark contribution xq(x) = xu(x) + xd(x) + xs(x) at Q 2 = 3 
GeV 2 (open circles) and Q 2 = 5 GeV 2 (full triangles) versus x. Data 
are from ref.[10] and solid line is the result of our fit. 

Fig. 4 The gluon distribution xG(x) at Q 2 = 5 GeV 2 versus x. Data are from 
ref.[ll] (area between dashed lines) and ref.[12] (area between small 
dashed lines) and solid line is the result of our calculation at Q 2 = 4 
GeV 2 . 

Fig. 5 xg[(x) at < Q 2 >= 10 GeV 2 versus x. Data are from ref.[5] and solid 
line is our prediction at Q 2 = 4 GeV 2 . 

Fig. 6 xg™(x) at < Q 2 >= 2 GeV 2 versus x. Data are from ref.[13] and solid 
line is our prediction at Q 2 = 4 GeV 2 . 
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